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In order to explain the present accelerating expansion of the universe and the re- 
lated dark energy with negative pressure, the scalar model and cosmological factor 
have been widely studied. Relatively, the spinor field is unreasonably overlooked. 
Considering the facts that all fermions are described by spinors and the uniform 
scalar field can hardly explain the galactic structure, so the dark spinors may be 
partially responsible for dark matter. In this paper, we give a detailed investigation 
on the state functions of spinors, such as the mass-energy density and equation of 
state, in the context of cosmology. Here the spinors are quantized and identified 
by nonlinear potentials, and their state functions are derived from standard statis- 
tical and variation principles. The results provide us some interesting consequences 
and new insights: (I) The negative pressure exists, and the equation of state can 
be designed by adjusting the nonlinear potential. (II) The initial singularity of the 
space-time can be removed by the nonlinear potential. (Ill) A short period of infla- 
tion exists and in this period the scale factor inflates much faster than the exponential 
function. 
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I. INTRODUCTION 

Since the discovery that the expansion of the universe is accelerating[l]-[6], the standard 
model of cosmology has shifted from a matter dominated and decelerating expansion picture 
to search for the dark energy and matter with repulsive gravity or modification of the general 
relativity. The present methods to analyze the observational data include direct measures of 
cosmic scales through Type la supernova luminosity distances, the angular distance scales of 
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baryon acoustic oscillation and cosmic microwave background density perturbations, as well 
as indirect probes such as the effect of cosmic expansion on the growth of matter density 
fluctuations. The basic building blocks of the universe is now believed to consist of 4% 
baryons, 20% dark matter, and 76% dark energy The cosmic accelerating expansion results 
in a profound mysteries in all of science, with deep connections to both astrophysics and 
particle physics. 

The simplest choice for such dark energy with negative pressure is the cosmological con- 
stant A, a term consists with Einstein's equations. This term acts like a fluid with equation 
of state(EOS) p\ = —p\. Although the ACDM model is basically consistent with all obser- 
vational data, the evaluated A seems incoincident with the local Newtonian approximation 
of the galactic gravity, and a lot of theoretical efforts was paid to explain its origin and 
apparent value [7] -[10]. 

The widely studied dynamical models to describe dark energy are mainly the scalar 
fields, such as quintessence, K-essence, tachyon, phantom and dilaton, which are uniformly 
distributed in the universe[10]-[19]. These models have different theoretical origin and con- 
sequences, and occasionally can give some explanation for the observational data and the 
origin of A. However these models, viewed objectively, suffer from several shortcomings 
such as lack of predictive power and the potentials without natural theoretical justification, 
etc. [11, 18, 20]. 

Noting the facts that all fermions are described by spinors and the uniform scalar field can 
hardly explain the galactic structure [21], so the dark spinors may be partially responsible 
for dark matter. In this paper we give a detailed discussion for their state functions, such 
as the mass density and equation of state, in the context of cosmology. Here the spinors 
are quantized and identified by nonlinear potentials, and their state functions are derived 
from standard statistical and variational method. This work is a further research of the 
preceding papers [22] -[25]. Some similar works were once done in[26]-[29], where the authors 
also considered the nonlinear spinor field as candidate of dark energy, but only one spinor 
field was taken into account. 

The nonlinear spinors provide us some interesting consequences and new insights: The 
negative pressure exists, and the EOS can be designed by adjusting the potential. The initial 
singularity of the universe is cancelled by the nonlinear potential. A short period of inflation 
without inflaton exists, in which the scale factor inflates much faster than the exponential 
function. The state functions are only the functions of scale factor, but independent of its 
derivatives. Combined with other models, the nonlinear spinors are hopeful to give more 
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natural explanation to the cosmic accelerating expansion and other observational data. 



II. FUNDAMENTAL CONCEPTS AND EQUATIONS 

At first, we introduce some standard concepts and equations in cosmology and field 
theory. In this paper, we adopt the conformal coordinate system for the universe, then the 
corresponding Friedmann-Robertson- Walker metric becomes[22] 

<7 M „ = a 2 (£)diag[l,-l,-S 2 (r),-S 2 (r)sin 2 6] , (2.1) 

where 

sin r if K = 1 , 
[ r if K = 0, (2.2) 
sinh r if K = — 1 . 

The Lagrangian of the gravitational field is given by 



where R is the scalar curvature, A the cosmological factor. For the metric (2.1), we have 
the nonzero components of the Einstein's tensor as follows 

<% = AK + k)> G\ = GI = GI = -\( 2 -^-^ + k), (2.4) 



a z \a z J a z \ a a z y 

where a' = j|a, S' = -£,S. Corresponding to metric (2.1), the energy-stress tensor certainly 
takes the following form 

T^ = diag[p,-F,-P,-F], (2.5) 

which is unnecessarily related to the perfect fluid model. In (2.5), p stands for the total 
energy of all matter, but P includes not only the usual pressure, but also all potentials of 
interaction [23, 24, 30]. 

By the Einstein equation, for the metric (2.1) we have the following two independent 
equations, 



= G° + A + 8nGp = ~^K 2 + Ka 2 ) + A + 8nGp, (2.6) 
= G£ + 4A + 8nG(p - 3P) = + Ka) + 4A + 8vrG(p - 3P), (2.7) 
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(2.6) and (2.7) lead to the energy-momentum conservation law 

1 d(pa 3 ) 
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3a 2 da 

If reckoning in the effect of A, we get the total energy-momentum tensor with vacuum energy 
r^|^ = diag[p to t,-^tot,-^tot,--Ptot] J Ptot = P+^^, Ptot = P--^—^- (2.9) 

07TG 07TG 

Correspondingly, the Einstein equation becomes + 87rGT^| tot = 0, and then the total 
EOS is defined by 

P tot 8ttGP-A 

w= = T- (2.10) 

p tot 8ttG> + A 

Now we construct the Lagrangian of the nonlinear spinors. Denote the Pauli matrices by 

*-^-{C:)-0o)-0-°0}' m) 

Define 4x4 Hermitian matrices as follows 

"{(:;)•(::)}-(-)•'•(-)• <»> 

We consider the following Lagrangian for the same kind spinors in flat spacetime[22, 23, 24] 

c = Efc <f>t ( a k id i* - ^7)0fc + F{%, Pk, ^ 13 ^ 

ik = <j>ti(pk, fa = ^iPA, &k = <Pia^<t>k, 

where we assign (pj, to the k-th dark spinor Sk, \i > is constant mass, which takes one 
value for the same kind particles. F is the nonlinear potential term. 

Considering that all the particle-like solutions of the nonlinear spinor equation have a 
mean diameter about 10 2 Compton wave lengthes, and all \(pk\ decay exponentially with 
respect to the distance from the center [23, 24]. So we can neglect the nonlinear interaction 
terms among spinors, but keep the self coupling terms. Then the nonlinear potential F 
should take the following decoupling form 

F = ^V fc , V k = V(%), (2.14) 

k 

where the self-coupling potential V is a differentiable and concave function. 
The generalized form of (2.13) in curved space-time is given by[31, 32] 



The variation of (2.15) with respect to §\ gives the dynamic equation for each S k 

It is easy to check the normalization condition holds for each spinor 

/ \<f> k \ 2 a 3 dtt = 1, (Vfc), dtt = S 2 sin 9drd9dip, (2.17) 
in 

where dfi is the comoving volume element independent of t. 

Coupling (2.3) with (2.15), we get the total Lagrangian of the universe with dark spinors 

L = J^(R ~ 2A) + E tttfidpfa) - fi% + V k ) y/g, (2.18) 

k 

where y/g = \J\ det(^)| = a A S 2 sin 9. Generally the variation of (2.18) with respect to g^ u 
gives the Einstein equation 

+ Ag^ + 8nGT^ u = 0, (2.19) 
where the energy momentum tensor T^ v is given by (only valid for diagonal metric [31]) 

T " U = E {WW + eW)&) + {Vi% - V k )gA . (2.20) 

k ^ ' 

By the dynamic equation (2.16), we get 

= ^(^ + 3^-4^). (2.21) 
By (2.5), the following relation holds 

P ~ 3^ = I* = 75 / E ^ + - 4Ffc ) ( 2 - 22 ) 

III. THE EQUATION OF STATE 

In this section, we derive the state functions p{a) and P(a) corresponding to the following 
nonlinear potential V, but the procedure is valid for more complex cases. 

v = ^ + ^ + - + ^-r» (3.1) 

2 3 n + 1 
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Denote the present time by to an d do — a (^o)- Making transformation 

_ 3 

*=-, - = (3.2) 
and considering that in the region of spinor particle cioS(r) = a S C^~) ~ > ^> we get 

V = Hi (2°)'? + !5 (!!>)> + . .. + ^ (2°)^ V, (3.3) 

2 \ a / 3 V a / n + 1 V a / 

V = „ 1 (5V I » ! (5 5 7 ! + ... + »„(^)> 1 (3.4) 
afidnip = ^-(a-wi (^j ^ - w 2 (^j 7 2 w n 7™^ 7^- (3-5) 

in which 7 = -0 + 7^. The normalization condition (2.17) becomes 

1 = y i0i 2 a 3 dfi = y i^i 2 a^fi = y mvz, (3.6) 

so ip is normalized in the new coordinate system. By (3.2)-(3.4), equation (2.21) becomes 

rpfj, Sr(~f a °Y> ~2f a o\ 6 . , ~n+i3jj7-l /a \ 3 ( n+1) \ , Q7 ^ 
T ^22{^{-) +^*(-J + - + «vy* "^(-J J- (3-7) 

(3.7) is a scalar equation, which is more convenient for local Lorentz transformation. We 
take (3.7) as an entry to derive state functions. 

Denote fij to be the static parameters defined in central coordinate system 

A*o = / LHkd 3 x, fij = / Wj% +1 d 3 x, (j = l,---,n). (3.8) 

Since a(t) is a slowly varying function, we can treat it as constant. Then all in (3.8) only 
very weakly depend on a, and they can be regarded as constants except for the case a — > 0. 
Assuming the peculiar speed of the k-ih spinor is v^, which depends on a(t). Solving the 
geodesic, we have [30, 33] 

Vh = — , ^ k , or \ 1 1 — vl = — ; a , (3.9) 

where b k is a constant determined by initial speed. By (3.9), we can calculate the mean 
value of (3.7) as 

- ^7 e + ^ er + • • • + (3» - iK ) 

A/"/ /a \ 3 x % /a \ 3 ™\ / 3a kT \ 

= * ("» + 2 '" (t ) + ' ' ' + (3n - 1K (t) j - ^Tlr) j • <" > 
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where Xk stands for the central coordinate of Sk, and dV the kinetic energy distribution, 
A/" = the comoving number density of the spinors, m is mean mass of all particles, a is a 
parameter determined by distributi 
In [33], we derived the relation Tia 



parameter determined by distribution function. For Maxwell distribution, we have a — |. 



omb am 



kT(a) = = — (v^T^ - a), (3.11) 

where h is a constant determined by the initial temperature. 

It should be pointed out that, (3.9) and (3.11) only exactly hold for mass points, but 
approximately hold for nonlinear spinors with a tiny mass and a large potential [30, 34]. This 
may be the underlying reason why the ordinary baryons in a galaxy can be automatically 
separated with dark matter and shift into a disc to develop the spiral structure[21]. However, 
in the case a ^> b and (vk — > 0,T — > 0), the following results hold to high precision. In the 
extreme case near a — > 0, only the trends of (3.8)-(3.11) are right. 

In order to get more information and intuition, we calculate the following simple case in 
detail. Assume 

q = A/>o, e = —al> 0, ^ = 0, (j > 2), (3.12) 

where g is the comoving density independent of t or a. Substituting (3.11) and (3.12) into 
(3.10), we get 




and then 



J a 3 T^da = -3ga + ga ^1 + ^-(V~c 
ge I 3cr yfa 2 + b 2 — 



b 2 



(3.14) 



a 2 \ 2 V /7#TF + 
On the other hand, the trace of Einstein equation (G^ + 4A + 8ivGT^) = gives 

8nG J a 3 f^da = 3(a /2 + Ka 2 ) - Aa 4 , 

= -a\Gl + A) = 8nGpa 4 . (3.15) 
Using (3.14), we solve the rigorous p(a) as 

p = h I a%T " da 



v l + ^(x^^-a)V^-^(l-^ V^~ a h (3-16) 
a 3 V 2a v ') a 4 a 6 I 2v^2^2 + a ' 



s 



where a is a constant probably related to the gravitational potential, the first term corre- 
sponds to the mass energy and kinetic energy of the spinors, and the second and third terms 
correspond to the total potential energy including the gravitational one. Substituting (3.16) 
into the energy-momentum conservation law (2.8), we get the rigorous P(a) as 



a 4 2aVa 2 + b 2 a 6 \ 2(Va 2 + b 2 + a)Va 2 + b 2 J ' l ' ' 

We can check that, (3.16) and (3.17) also satisfy the equations (2.7) and (2.22), so they are 
the exact cosmic EOS for the spinors. The constants (g,a,e,b,a ~ |) can be treated as 
macroscopic parameters to be determined by empirical data. However, it should be keep in 
mind that, the above results have only qualitative significance when a — > 0. 

The constant a is equal to — |Co in [22], where the explanation for its physical origin and 
meaning is superfluous and inexact. By (3.16) and (3.17), we find 

P a 

► >0, (a->oo), (3.18) 

p a 

which proves the infinitesimal value of the dimensionless pressure for a spinor in flat space- 
time as worked out in [24]. 



IV. GLOBAL PROPERTY AND ASYMPTOTIC BEHAVIOR 

Substituting the above EOS into (2.6), we get a closed dynamic equation for scale factor 
a(t), 

3ab 2 e ( 3a 



a ' 2 = - Ka 2 + + 2 R 
= -$(a), 



3a k 1 £ (a 

2(V?+F+a) a 2 V 2 ^ v ' 



(4.1) 



where R is the mean scale with length dimension. 

r = ^, y^T£^a_ (42) 

Of course, for practical resolution to the cosmological model, one should introduce the 
mass-energy density of other fields such as photons to the total potential $(a). If taking 
a(t) as generalized coordinate of a mass point, then (4.1) have the form of Hamiltonian, and 
$(a) is equivalent to the potential. The qualitative analysis for its solution is routine and 
convenient [19, 35]. 

Since e > 0, by (4.1) we find that, if the nonlinear potential of the spinor exists, the initial 
singularity is absent due to $ — » +oo as a — > +0, namely, the phase trajectory of (a', a) 
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can not reach a — > 0. In the cases of ordinary fluid model and linearized Dirac equation, we 
have actually assumed e = 0, and then the fact is concealed. 

For the behavior of solution of (4.1) and the destiny of the universe, the roots of $(a) = 
have decisive significance [3 5]. Since e > 0, we have at least one positive root. Denote the 
least positive root by Co, the second positive root by Ci(Ci < +°°)) where Ci = +oo means 
the unbound trajectory and bouncing universe, but Ci < +oo corresponds to the cyclic 
universe. Then we only need to discuss the solution in the connected region Co < a < Ci- 

Near the root C = Co or C = Ci < 00 > we have the following Taylor expansion 

a' 2 = A (a -Q + A^a- () 2 + 0((a - C) 3 ). (4.3) 
In contrast it with (4.1), we get the parameters expressed by C as 

e ( 3a ^/e + ¥-C \ 1 oh 2 C -3KC 2 + C 4 A 

3C 2 ^ 2 A /c2T^ + c; + 2c+ v ^T^ + 3 + 18R 




C 4 \ \ V(C 2 + b 2 ) 3 J J 2v/(C 2 + 6 2 ) 3 
The solution near the roots equals to 



1 



a = ( + -A t 2 + O(t A ). (4.7) 



In the case C - ► Ci < 00 ? we should have Aq < 0. So for the cyclic universe, by (4.5) and 
Aq < 0, we have the estimation K > |ACi + ^ > 0, which implies K — 1, and then 



3 Ci - i2) 2 

A < To < —77. 4.8 

2d 3 ~ 9i? 2 



In the case C ~~ ^ Co ~ +0, we should have t4 > 0, omitting the higher order terms, then 



we have 



a > o, if e < A - ^ vg+g-c»\ C 5 + c \ c (4 . 9) 

If a > 0, the trends of all functions are similar to those displayed in [25]. 
Denote r = f Q a(t)dt to be the proper time, by (4.3) we get 

d 2 a A R . 

J-2 =2/5-73-°°, (0,-0). 

aT T =0 Zt >0 SO 
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So the inflation indeed exists for any nonsingular cosmological model, and the scale factor 
would initially expand much faster than the exponential function. The temperature param- 
eter b have little influence on the global behavior, because terms including it always have a 
large part, but it has important influences on the value of parameters (Co, 

Although the above w(a) for the nonlinear spinors is negative near £ , the trends is 
different from that of the fitted EOS previously derived from the observational data. The 
former is an increasing function of a in the case of a > 0, but the latter seems prefer to 
a decreasing function. However, according to the WMAP five-year data[6], if dark energy 
perturbation is included, the constraints on the time evolving equation of state of dark 
energy is consistent with an increasing u>(a)[36]. 

On the other hand, similar to the quintessence model, by choosing suitable potential, we 
can almost get any wanted EOS [18]. Calculation shows that, if the potential of the spinors 
including terms like W-i(C 2 + 7 2 ) - ^, the corresponding EOS is decreasing function similar 
to the fitted functions. But few people would believe the Nature works in such manner. 

V. DISCUSSION AND CONCLUSION 

According to the above calculation, we get some important consequences and insights. 

1. By (3.17), for the nonlinear spinors, the negative pressure definitely exists, and the 
equation of state can be designed by adjusting the potential. 

2. By (4.1), we find the initial singularity of the universe is cancelled by the nonlinear 
potential of the spinors due to <3>(a) < when a — > 0. 

3. A short period of inflation of the universe exists and in which the scale factor inflates 
much faster than the exponential function[25]. Such inflation implies the universe only 
spends a very short time to cross the turning point. 

4. Although the dynamical equation of the spinors is more complex than that of the 
scalar fields, the state functions of the spinors are even relatively simpler due to the 
normalizing condition (3.6). These functions are only the functions of scale factor, but 
independent of its derivatives. 

5. By (3.16), we find that, different from the point particle and fluid model, mass-energy 
density of the spinors is not a monotone decreasing function of a. Spinors have strong 
interaction with the heavily curved space-time. 
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